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New millimeter-wavelength emission lines are reported for the Cloverleaf 
(H1413+117), a gravitationally lensed quasar at redshift 2.56. High signal-to-noise 



£ — \ ratio spectra of four lines in CO (J = 3 — > 2, 4^3, 5 — > 4, and 7 — > 6) have been 

obtained and are modelled using an escape probability formalism. The line brightness 

^ ' temperatures are flat or rising between CO(3 — » 2) and CO(4 — » 3) and then drop to 

the CO(5 — ► 4) and CO(7 — ► 6) lines; this falloff suggests that the optical depths in the 
CO lines are modest (t4_3 < 3), and that the gas is relatively warm (T ^ 100 K) and 
dense (n Ha ~ 3 x 10 3 cm -3 ). 

The neutral carbon fine structure line CI( 3 Pi — ► 3 i 3 o) was detected with moderate 
SNR on two separate occasions, and appears secure. The HCN(4 — > 3) line is somewhat 
more tentative, as it was observed only once and detected with moderate SNR. 
However, its brightness ratio relative to the CO(4 -» 3) line, L'(HCN)/L'(CO) « 0.25, 
is very similar to the ratio seen in luminous infrared galaxies. This probably reflects 
the gas fraction maintained at high densities by radiation pressure from the active 
nucleus; alternatively, the HCN transition may be excited through absorption of 
infrared radiation at 14^m. 

The detectability of the Cloverleaf in molecular and neutral atomic transitions is 
largely due to high emissivity of the gas and magnification from gravitational lensing, 
rather than an extremely large mass of gas. 

Subject headings: quasars: individual (H1413+117) — ISM: molecules 
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1. Introduction 

There are few observational windows on galaxies at high redshift, and each new one that 
opens up is invaluable. It is now possible to directly study interstellar gas in significant quantity 
in three systems at z > 2: IRAS F10214+4724 (z = 2.28), the Cloverleaf quasar H1413+117 
(z = 2.56), and the quasar BR1202-0725 (z = 4.69). F10214+4724 was first identified as a 
luminous infrared galaxy (Rowan- Robinson et al. 1991), and has been detected in three transitions 
of carbon monoxide (Brown & Vanden Bout 1992a; Solomon, Downes, & Radford 1992a). It is 
now known both that this object is gravitationally lensed (e.g. Broadhurst and Lehar 1995), and 
that it harbors an obscured quasar nucleus seen in polarized (reflected) light (Goodrich et al. 
1996). BR1202-0725, the third highest redshift quasar currently known, has recently been detected 
in CO by Ohta et al. (1996) and Omont et al. (1996). The Cloverleaf quasar, the subject of this 
paper, is a QSO with classical broad emission lines plus broad absorption troughs in the optical, 
with an image that is lensed into four spots separated by about 1" . 

Based on its strong submillimeter continuum emission, indicative of a large mass of cool dust 
(Barvainis, Antonucci, &; Coleman 1992), the Cloverleaf presented a good candidate for the study 
of molecular gas at high redshift. Indeed, CO(3 — > 2) line emission, redshifted to 97 GHz, was 
searched for and readily detected (Barvainis et al. 1994; Wilner, Zhao & Ho 1995). The close 
similarities in both infrared/submillimeter continuum and CO emission properties between the 
Cloverleaf and F10214+4724, and the fact that both have powerful active nuclei, have given rise 



to models uniting these two objects through orientation/obscuration effects (e.g., Barvainis et al 



1995; Granato, Danese, & Franceschini 1996), and have bolstered evolutionary scenarios linking 



luminous infrared galaxies and quasars in general (e.g., Sanders et al. 1988; Wills et al. 1992). 

To provide more information on the state of the interstellar medium in the host galaxy 
of the Cloverleaf, we have obtained detections in three new transitions of CO, and have made 
measurements in two transitions of neutral carbon and one of HCN. We present these results, 
together with models for the emitting gas, below. 



2. Observations and Results 

The observations were carried out using the IRAM 30m telescope on Pico Veleta, Spain, 
during observing runs in 1994 June (two nights) and 1995 May (seven nights), and with the 
Plateau de Bure Interferometer (PdBI), near Grenoble, France, during 1996. At the 30m the 
weather was mainly clear during both runs, but occasional high winds during the 1995 run 
resulted in some downtime. Spectral data were obtained in the 3mm, 2mm, and 1mm bands 
simultaneously, using SIS mixer receivers and filterbank or acousto-optical backends of 512 MHz 
bandwidth. System temperatures derived using the chopper wheel method were typically 250 K, 
200 K, and 400 K in the 3mm, 2mm, and 1mm bands respectively. Double position switching was 
employed to provide flat baselines, with subreflector nutation at a rate of 0.5 Hz. 



Clear detections were obtained in four CO transitions: J = 3 — ► 2, 4 — ► 3, 5 — > 4, and 7^6. 
Peak antenna temperatures were in the range 4-8 mK. 13 CO(3 — ► 2) was searched for but not 
detected. A similar line flux and velocity width for the J = 3 — ► 2 line was obtained by Wilner et 
al. (1995). The C0(7 — > 6) single-dish measurement presents a problem, because of the limited 
available bandwidth of 512 MHz. Corresponding to 680 km s _1 at 226 GHz, this is not wide 
enough to provide a zero-level baseline - even with subreflector nutation the zero level in the 
raw spectrum cannot be trusted. We will therefore instead use here the PdBI observations of 
CO (7 — ► 6) (for details, see Alloin et al. 1997). These also suffer from the same limited bandwidth, 
but for interferometer observations the zero level should be accurate. To derive the total intensity 
(Jy km s _1 ) of the line for comparison with the other CO transitions it is necessary to fit the line 
with a gaussian, but the line width is poorly constrained since the line wings are cut off at the 
band edges. Fixing the linewidth to the mean of the other three CO transitions, AVfwhm = 376 
km s _1 , we obtain IpdBi = 43.7 ± 2.2 Jy km s _1 (the error is statistical only, and does not include 
10% calibration uncertainty). 

The spectra are shown in Figure 1, in units of flux density. Table 1 gives line parameters 
and upper limits as derived from Gaussian profile fits to the 30m spectra [except for CO (7 — ► 6), 
where the line parameters are estimated as described above]. The line luminosities are given 
in two ways: L, in units of solar luminosities, and V in units of K km s" 1 pc 2 (see Solomon, 
Downes, & Radford 1992a for definitions and formulae for computation). 1/ is proportional to line 
brightness (Rayleigh-Jeans) temperature integrated over the area of the source, and assuming that 
the CO lines arise in the same general volume this quantity can be used to form line brightness 
temperature ratios to constrain physical conditions in the gas (see below). 

Observations were also made of the neutral carbon fine structure lines CI( 3 Pi — ► 3 Pq) and 
CI( 3 P 2 -► 3 -Pi) (rest frequencies 492.161 GHz and 809.343 GHz, respectively). A detection in the 
former made in the 1994 run repeated with reasonable consistency in 1995, and the average of all 
good scans from both years is shown in Figure 1. Confidence in this detection is enhanced by the 
similarity in redshift and line shape to the CO lines (Table 1). Also, the line repeated after a shift 
in observing frequency which placed it in a different part of the filterbank. In CI( 3 i-2 ~~ ► S Pi), only 
an upper limit was obtained (Table 1). 

Finally, we observed HCN(4 — ► 3) during the 1995 run. We report a tentative detection of a 
line with peak flux density of 5.5 mJy (Figure 1). The total integration time for this observation 
was about 30 hours. Although this line is near the limit of detectability for the 30m telescope, as 
in the case of CI( 3 Pi — ► 3 Pq) we are encouraged by the similarity of the line parameters to those 
of CO. 

Some comments on the uncertainties in the observed line ratios are in order. In addition 
to the formal errors from noise in the measurements, it is necessary to include the systematic 
errors. Conversion from antenna temperature in Kelvins to flux density in Janskys was done 
using standard calibration factors for the 30m telescope derived from observations of planets and 



other standard sources. The errors in these factors are estimated to be ±10% for the 3 — ► 2 
line and ±15% for the 4 — > 3 and 5 — ► 4 lines (C. Kramer, private communication). For the 
PdBI observations of CO(7 — ► 6) the calibration uncertainty is estimated at 10%. Including the 
systematic and statistical errors, the CO line Rayleigh-Jeans brightness temperature ratios relative 
to the strongest line, CO(4 -> 3), are then: (3-2)/(4-3) = 0.83±0.16, (5-4)/(4-3) = 0.73±0.16; 
(7 — 6)/(4 — 3) = 0.68 ± 0.13. These ratios are computed from the V values given in Table 1. 



3. Interpretation 

3.1. The CO lines 

The CO lines are characterized by a constant or increasing brightness temperature in going 
from J = 3^2toJ = 4^3, followed by a decrease in the higher-J transitions. We cannot 
rule out the possibility that the brightness temperatures are the same for all transitions, since the 
line ratios differ from predicted thermal R-J values by only 1.5 — 2er in each case. In the case of 
thermal line ratios we would be able to say little about the column density or temperature of the 
gas. However, below we will assume that our measured values represent the true line brightnesses, 
and use them to model physical conditions in the gas. This exercise suggests that the line optical 
depths are modest, with r ~ a few or less. 

In modeling the CO line ratios, we calculate the excitation of the CO molecules using a 
non-LTE scheme which employs mean escape probabilities to include the effects of radiative 
trapping on the excitation of the CO (cf. [Aalto et al. 1991 ); the results depend on the density n 



of H2, the gas kinetic temperature T, and the CO column density per unit linewidth Nqq,/ AV . 
We note in passing that in non-LTE models using the essentially equivalent LVG approximation, 
the ratio of CO abundance to velocity gradient - which fixes the optical depths - is often set to a 
value appropriate to Galactic molecular clouds (e.g., Solomon, Downes, & Radford 1992a), which 
is not a priori justifiable when the physical conditions in the gas may be radically different. 

In Figure 2 we show the allowed ranges in n and Nqq,/ AV for three temperatures, T = 60, 
100, and 300 K; the allowed ranges are the hatched regions within the curves. A temperature 
of 60 K is the minimum acceptable temperature. We have also imposed the constraint that the 
brightness temperature of the CO(7 — > 6) line must be at least 40 K.f] This constraint can in 
principle eliminate a number of models which satisfy the line ratios. For the high temperature 



2 Using the observed lensed CO source size in the CO(7 — > 6) line (Alloin et al. 1997) indicates a source radius 
r < 0A8h~ 1 rriQ 1 kpc (assuming go = 1/2), where mc ~ 0.2m is the magnification of CO in image C of the quad 
(m is the total magnification for all four images). The exact relation between mc and the estimate of the true 
source radius r is model dependent - here we have assumed that lensing stretches the source predominantly in one 
dimension; see Burke et al (1991) for an illustration of distortions expected in a quad lens. Given these assumptions, 
the minimum brightness temperature in the CO(7 — » 6) line (that is, for unit areal filling factor) is Tr(7 — 6) ^, 40mc 
K, independent of h. 



models, the allowed range in density is a factor of 3 — 10, while for T = 60 K any density greater 
than n ~ 8 x 10 4 cm -3 is allowed. The acceptable range of Nqo/ AV increases from about a factor 
of two (T = 60 K) to a factor of seven (T = 300 K). 

At T = 60 K, the line ratio constraints can be satisfied for thermalized level populations 
(which is why there is no upper limit to the acceptable density), provided that the optical depth 
is not too large (t3_2 £ 2): otherwise the line ratios would simply be the ratio of Rayleigh- Jeans 
brightness temperatures for T = 60 K, which is too large compared to the observed line ratios.(j 
Thus the ratio of the optical depth terms (1 — e~ T ) must be less than one, which requires that the 
optical depths are not much greater than unity and that tj^q is less than t-^2- The sharp cutoff 
in Nqo/ 'AV at the high end is the result of this optical depth restriction; the boundary is nearly 
vertical because, except at the lowest densities, the rotational levels are thermalized through J = 7 
for all the densities. 

At high temperatures the line ratios are matched by somewhat different physics. Thermalized 
levels are unlikely, as the high— J lines would then be too strong relative to J = 4 — > 3 no matter 
what the optical depths. The observed fall-off to the higher J lines can still be matched in this 
case, provided that the level populations are subthermal; in particular, the excitation temperatures 
T ex of the J = 5 and J = 7 levels must be well below that of J = 4. This restricts the density from 
being too large, and in fact there is only a narrow range (at a fixed temperature) at which the 
line ratio constraints are satisfied. In addition, the optical depths are also restricted from being 
too large, as otherwise radiative trapping would inhibit the required decrease in the excitation 
temperature. It is really the drop in T ex that is crucial, as the optical depths may actually increase 
from the 3 — > 2 line to the 7 — > 6 line; the drop-off in line flux is then due to the fact that kT ex /Ej 
is less than ~ unity for the J = 7 level, while T ex (J = 3) is non-negligibly larger than T ex (J = 7). 
There is a systematic trend in the line ratios; the model values approach the maximum allowed 
observational values with increasing iVco/AV. 

The maximum optical depth in the CO (4 — > 3) line is generally ~ 3. The relatively low 
optical depths are mandated by the low value of the (5 — 4)/(4 — 3) and (7 — 6)/(4 — 3) ratio, 
as outlined above. This result is unusual in comparison to Galactic molecular clouds, where the 
optical depth in the CO J = 1 — > line (the most frequently observed transition) is generally 
inferred to be large. However, modest optical depths in the CO rotational transitions appear to 
be common in starburst galactic nuclei flAalto et al. 1995 ), for the same reason as in our models 



of the CO emission from the Cloverleaf: the elevated gas temperatures (T ^ 60 K) distribute the 
molecules over a much larger number of rotational levels. In addition, since the Cloverleaf is at 
a redshift z = 2.558, it is possible that the gas-phase metallicity is considerably below the Solar 
neighborhood value, which will affect the CO abundance and emission (e.g. Maloney & Wolfire 
1996). A clear prediction of our single-component CO models is that the CO(l — » 0) line should 



3 For example, the R-J brightness temperature ratio for optically thick and thermalized lines at 60 K is 0.89 for 
(7-6)/(4-3), compared with the observed value of 0.68 ± 0.13. At 300 K, the thermal R-J ratio is 0.98. 



be weak compared to the CO(3 — ► 2) line, with a brightness temperature about one-third that of 
the 3^2 line; this is a consequence of the increase in optical depth with J up to J ~ 4 or 5. 



3.2. The CI lines 

Our upper limit on the intensity of the CI 809 GHz line gives an upper limit to the 
CI( 3 i-2 — * 3 -Pi)/CI( 3 Pi — > 3 i"o) ratio of ~ 3.9 (here using luminosities L expressed in L@). We 
have calculated the expected ratio of these lines as a function of density, T, and Nc/AV, as for the 
CO lines. The upper limit to the 809/492 ratio is only significant for gas temperatures T £ 100 
K, where it restricts the gas densities to be n ;$, 10 4 cm -3 . Thus in principle, the CI emission can 
arise from the same regions which produce the CO emission, as delineated in Figure 2. 

The observed ratio of luminosities L in the CI 492 GHz and CO (3 — > 2) lines is approximately 
0.5 (Table 1). In order to match this ratio, the column density per unit linewidth of neutral carbon 

Therefore, if the CO and CI emission arise from 



-i\-i 



must be N C /AV > 3 x 10 17 cm" 2 (km s 
the same regions, the neutral carbon abundance is comparable to or a few times larger than that 
of CO. This unusually high ratio of C to CO is a reflection of the rather low ratio of Nqo/AV 
required to match to CO line ratios. 

We can compare the CI/CO ratio for the Cloverleaf with the COBE results for the Milky Way 
(Wright et al. 1991), and extragalactic detections of the 492 GHz line in IC 342 by Biittgcnbach 
et al. (1992)| , and in M82 by [Schilke et al. (1993)| and |White et al. (1994)fl . Relative to the CO 
J = 2 -► 1 line, the flux in the 492 GHz CI line is 1.4 ± 0.4 (IC 342) and 2.3 ± 0.6 (Milky Way). 
If the CO lines are optically thick, the ratio of the CI line to the CO transitions will scale as v~ 3 
in the high-temperature limit, so that relative to the CO (3 — > 2) line these flux ratios are 0.4 (IC 



342) and 0.7 (Milky Way), similar to the value of 0.5 seen in the Cloverleaf. From Schilke et al.| 
|(1993)| and [White et al. (1994)| , the ratio in M82 is 0.7. [Biittgenbach et al. (1992)1 conclude that 
the ratio seen in IC 342 and the Milky Way is consistent with the bulk of the emission arising in 
photodissociation regions (PDRs), although [Schilke et al. (1993)1 and [White et al. (1994)1 argue 



that the ratio is too large compared to standard PDR models; the CI intensity can be enhanced 
if the PDRs are clumpy or the cosmic-ray flux is very high. Although the CI and CO emission in 
the Cloverleaf probably do not arise in PDRs (see §3.3 and §4) the similarity in the line ratios is 
probably a reflection of a similarity in the physical conditions: relatively high temperatures and 
densities, optically thin CI emission, and (marginally) optically thick CO emission. 



4 Brown & Vanden Bout (1992b) have reported detection of CI in F 102 14+4724, but at a substantially different 
velocity from CO (3 —* 2) . The emission in the two species is therefore unlikely to be cospatial, so we do not cite a 
CI/CO ratio for this object here. 



3.3. The HCN line 

Our tentative detection of HCN(4 — > 3) indicates a luminosity V (in K km s _1 pc 2 ) relative 
to the C0(4 -> 3) line of about 0.25, which is similar to the HCN(1 -> 0)/CO(1 -> 0) ratio seen 
in ultraluminous galaxies ( Solomon, Downes, Sz Radford 1992b[) , Except for T ~ 60 K, the lowest 



allowable value, the maximum allowed densities in the CO-emitting region (n ~ 10 5 cm -3 ) are far 
below the effective critical density required to thermalize the HCN(4 — > 3) line [n ~ 10 cm -3 ; even 
for the HCN(1 — > 0) line, the critical density is n ~ 10 5 cm -3 at these temperatures]. In this case 
the ratio of HCN and CO luminosities is essentially determined by the sub-thermal excitation of 
the HCN. (This is also true for models with T = 60 K unless the density is extremely high.) This 
is not a very satisfactory explanation, however, as it requires fine-tuning of the gas parameters 
to match the observed ratio, and does not explain why a ratio of order 0.2 is characteristic of 
ultraluminous galaxies. A similar constraint arises from the failure of Wilner et al. (1995) to 
detect the HCO + (4 — 3) line, implying that most of the gas is not at densities above the effective 
critical density for this transition, n cr ~ 2 x 10 6 cm -3 . 

It is more probable that this ratio reflects the fraction of the CO-emitting region which has a 
density high enough to collisionally excite the HCN line. Solomon, Downes, & Radford (1992b) 
argue that only star- forming regions in molecular clouds have such high gas densities, and that 
therefore ultraluminous galaxies are powered by star formation. However, this argument overlooks 
the fact that in galaxies with powerful active nuclei, such as the Cloverleaf and F10214+4724, 
there is a large minimum gas pressure which is enforced by the radiation pressure from the AGN.[] 

Assuming that the incident flux from the AGN is absorbed in a region which is thin compared 
to the thickness of a cloud, then in equilibrium the cloud pressure must at minimum equal the 
absorbed radiation pressure. In terms of P = nT, this constraint can be written 

UP -^44 ( jabs ~\ _-3 

tL V 0.1 



P min = 2 x 10 1U ^ ^ cm- J K (1) 



pc 

where the bolometric luminosity of the AGN is Lbol = 10 44 L44 erg s _1 , the distance of the cloud 
from the AGN is r pc parsecs, and the fraction of the incident luminosity which is absorbed is / a bs- 
In terms of the gas density rather than P, equation (1) is 

n min = 2 x 10 8 ^ (&f) T, 1 cm" 3 (2) 

where the gas temperature T = IOOT2 K. The bolometric luminosity of the Cloverleaf is 
-^boi ~ W^m^ 1 ^ 2 Lq, where m is the magnification factor. Modelling of the Cloverleaf images 



5 Here we assume the gas communicates the radiation pressure through much of the molecular column in directly 
illuminated clouds. This wouldn't necessarily be true with a highly imhomogenous medium with an effective covering 
factor far above unity. 



suggests that m ~ 10 (Kayser et al. 1990) but this is uncertain. Equation (2) predicts that the 



density will exceed 10 cm 3 [enough to collisionally excite the HCN(1 — ► 0) line] within a radius 

r = 8.8xl0 2 (^)" 1/a 3 5 - 1/8 (^) 1/2 pc (3) 

For n ^j 10 7 cm -3 , needed to excite the HCN(4 — ► 3) line, the numerical coefficient is smaller by a 
factor of ten. The apparent line fluxes for transitions which require very high densities to excite 
are likely to be significantly affected by differential magnification, as they will arise preferentially 
at small radii. 

An alternative possibility which may be generically important in infrared-luminous galaxies 
- either starbursts or AGNs - is that the HCN levels are not collisionally excited, but instead 



are populated through absorption of 14^m photons in the degenerate bending mode (e.g., Aalto 



et al. 1994 ). For optically thin HCN emission (no radiative trapping), the criterion for radiative 
pumping of the HCN(4 — ► 3) line to be significant is that the dust temperature fulfills the condition 

T d > 1025[ln(487VK + l)]^ 1 K (4) 

where the continuum has been assumed to be a Planck function multiplied by a dilution/optical 
depth factor W. For W = 1, equation (4) requires T d > 166 K; for W = 0.1, T d > 262 K. 
Substantial amounts of dust at these temperatures are undoubtedly present in the Cloverleaf 
( Barvalnis et al. 1995|) . We note that these two mechanisms predict rather different size scales 



for the HCN emission; from equation (3), radiation pressure will only keep the gas density high 
enough to collisionally excite the HCN(4 — > 3) line to a radius of 100 pc, whereas for a bolometric 
luminosity of 10 13 Lq, the dust temperature will be above 250 K to more than 300 pc distance. We 
might also expect that in the former case, the HCN and CO line profiles could differ significantly; 
testing this will require much higher signal-to-noise observations than are presently available. 



4. Discussion and Conclusions 

The presence of a luminous active nucleus appears to have profoundly influenced the state of 
the interstellar medium in the Cloverleaf. The high line luminosities, the high ratio of CI to CO 
line intensities, and the prominent HCN emission may be the result of the impact of the radiation 
from the AGN on the surrounding molecular clouds. 

The mass of gas in the Cloverleaf can be estimated from both the CO line observations and 
the CI 492 GHz line. The total molecular gas mass can be derived from an observed CO line 
luminosity and our detailed model results as 



10- 17 (iVco/Ay ) 
(yco/3 x 10- 5 ^ 



M H2 (CO) = L co (M cl /L cl ) « 2lbL co - ) 'J W 1 M Q (5) 



where Nqo/AV is in cm 2 km 1 s, v is the rest frequency of the transition, I v is the intensity 
at line center emerging from a cloud, and yco is the CO abundance by number relative to total 



hydrogen; M c \/L c \ is the mass-to-CO luminosity ratio of the clouds from our models. Using the 
observed luminosity of the CO(3 — ► 2) line of Lqo{2> — 2) = 1.0 x 10 8 m _1 /i -2 Lq, we find a 
molecular gas mass Mjj 2 ~ 2 x 10 10 ?ti _1 /i _2 Mq, essentially independent of the choice of model, 
(except for the highest column density, high-temperature models, where it can be as large as 
Mh 2 ~6x 10 10 ?7i _1 /i _2 Mq). This is an order of magnitude smaller than the value originally 



suggested by Barvainis et al. (1994) , due to the considerably higher emissivity of warm CO of 



modest optical depth. 

To estimate the gas mass from the CI emission, we assume that the 492 GHz line is optically 
thin. Over the allowed range of physical conditions (see §3.2) approximately 40% of the carbon 
atoms are in the J = 1 level of the ground 3 P state. Using an Einstein ^4-coefficient of 7.9 x 10~ 8 
s _1 ( Nussbaumer 1971] ), we find that the mass of hydrogen associated with the CI is then 



M ffl (CI) ~ 1400 ,fT n _ 5 , M Q (6) 

(yc/3 x 10 d ) 

where the line luminosity L492 is in Lq and yc is the carbon abundance by number. (In both 
equations [5] and [6] we have normalized to a carbon abundance of 1/10 solar; equation [5] assumes 
all carbon is in CO.) The observed line luminosity L492 = 5.8 x 10 7 m _1 /i -2 Lq, and so we obtain 
a gas mass Mm (CI) «8x 10 10 m _1 /i -2 Mq. 

If irradiation by the central source is important in pressurizing the clouds, there will be 
associated dynamical effects, as the radiative momentum imparted to the clouds will provide an 
outward acceleration. Scoville et al. (1995) suggested that support of the molecular gas against its 
own self-gravity by radiation pressure could be important in F10214+4724, in order to reconcile 
the fact that their derived molecular mass exceeds the dynamical mass (estimated from the 
measured CO size and linewidth) by an order of magnitude. We estimate Md yn for the Cloverleaf 
from the upper limit to the CO source size r (see footnote 2) and the observed linewidths: 

M dyn « ^?™ » 1.6 x lO 11 ™- 1 /*- 1 Mq, (7) 

G sin i 

where the inclination angle i has been taken to be 45°. Comparison with the above expressions 
for the masses derived from CO and CI shows that the gas mass is in both cases smaller than the 
dynamical mass (assuming, of course, that r is not much smaller than the current upper limit, 
and that the carbon abundance is not much less than 1/10 solar). Although radiation pressure 
support could be important in the Cloverleaf, it is not at present necessary to appeal to it. 

The gas mass derived from the CI observations is ~ 1.3 — 4 times larger than the molecular gas 
mass estimated from CO; either the hydrogen associated with the atomic carbon dominates the 
gas mass (plausibly a result of the high X-ray ionization and heating rates within a few kpc of the 



active nucleus: Maloney, Hollenbach & Tielens 1996), or else the CI emission is being differentially 



magnified by the lensing relative to the CO. In either case, the gas mass is substantial, but not 
extreme, especially as the total magnification m may be a factor of ten. The detectability of the 
Cloverleaf in CO and CI would appear to be the result not of an enormous mass of gas, but rather 



10 



the effects of gravitational lensing and the rather unusual physical conditions (high emissivities) 
in the gas, in consequence of the proximity of a luminous active nucleus. 

We are grateful to the referee, David Wilner, for insightful comments, and to the IRAM staff 
for excellent support during the observations. Stephane Guilloteau helped produce the CO(7 — > 6) 
spectrum from PdBI. PRM acknowledges support through the NASA Long Term Astrophysics 
Program under grant NAGW-4454. 
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FIGURE CAPTIONS 

Fig. 1. — CO, CI, and HCN lines from the Cloverleaf quasar, shown with Gaussian fits to the 
line profiles. See Table 1 for line parameters derived from the fits. The CO (7 — > 6) parameters 
have been estimated by fixing the gaussian-fitted line width to the mean of the other CO lines, as 
described in §2 of the text. 

Fig. 2. — Allowed regions of CO column density per unit linewidth versus molecular hydrogen 
density for three temperatures, based on non-LTE escape probability modelling of the CO line 
ratios. See text for discussion. 
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Table 1. Millimeter Line Observations of the Clover leaf. a 



Line 




^obs b 


1 A 


s v 


AVpWHM 


V c 


/ 


L/10 8 


V /10 10 






(GHz) 


(mK) 


(mjy) 


( km s _1 ) 


( km s~ 1 ) 


(Jy km s -1 ) 


(Lq) 


(K km s _1 pc 2 ) 


CO(3 


-2) 


97.199 


4.2(0.3) 


25.8(1.6) 


362(23) 


-29(12) 


9.9(0.6) 


1.8(0.1) 


13.4(0.8) 


13 CO(3 - 


-2) 


92.924 


< 1.1 


< 6.7 






< 2.4 


< 0.4 


< 3.6 


CO(4- 


-3) 


129.576 


7.8(0.3) 


52.7(2.0) 


375(16) 


+12(7) 


21.1(0.8) 


5.1(0.2) 


16.1(0.6) 


CO(5- 


-4) 


161.964 


7.2(0.4) 


55.6(3.2) 


398(25) 


+5(12) 


24.0(1.4) 


7.2(0.4) 


11.7(0.7) 


CO(7- 


-6) 


226.721 




109(5) 




-5(13) 


47.3(2.2) 


18.3(0.9) 


10.9(0.5) 


CI( 3 Pi -> 


3 Po) 


138.351 


1.1(0.1) 


7.7(0.8) 


430(46) 


-47(24) 


3.6(0.4) 


0.9(0.1) 


2.4(0.3) 


CI( 3 P 2 -> 


3 Pi) 


227.510 


< 2.2 


< 22.2 






< 8.5 


< 3.6 


< 2.1 


HCN(4 - 


-3) 


99.631 


0.9(0.2) 


5.5(1.1) 


436(103) 


+49(40) 


2.6(0.5) 


0.5(0.1) 


3.4(0.6) 



a Line parameters derived from Gaussian profile fits. Quoted uncertainties are statistical only. For CO (7 — > 6), the line 
strength was estimated from measurements taken at the Plateau de Bure Interferometer assuming a line FWHM of 376 km 
s — 1 , as described in §2 of the text. Upper limits are 3<r, and were derived from Gaussian fits with the linewidth and redshift 
fixed to the mean values for the detected lines. Line luminosities assume Hq = 75 km s~ x Mpc -1 and qo = 0.5, and are not 
corrected for gravitational magnification; see text. 

b Observed frequency of line centroid. 

c Vclocities measured relative to 2 = 2.5579, which is the weighted average redshift of the four detected CO lines. 
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